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ABSTRACT 

This report reviews theoretical studies performed by the Aeroballis- 
tics Division with regard to the ablation-type heat protection of vehicles 
traversing an atmosphere. These studies rest on numerical solution 
methods of the pertinent partial differential equations representing 
balances of mass, momentum, and heat. Provided sufficiently small steps 
are used, the numerical calculation methods yield exact solutions for 
homogeneous and noncharring materials. This report presents major con- 
clusions from these ablation studies; the details appear in twenty- three 
papers published by Aeroballistics Division between 1959 and 1963. 
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STUDIES ON ABLATION OF OBJECTS TRAVERSING AN ATMOSPHERE 
By Ernst W. Adams, John D. Warmbrod, and Benton K. Berry 

SUMMARY 

This report reviews theoretical studies performed by the Aeroballis- 
tics Division with regard to the ablation-type heat, protection of vehicles 
traversing an atmosphere. These studies rest on numerical solution 
methods of the pertinent partial differential equations representing 
balances of mass, momentum, and heat. Provided sufficiently small steps 
are used, the numerical calculation methods yield exact solutions for 
homogeneous and noncharring materials. This report presents major con- 
clusions from these ablation studies; the details appear in twenty- three 
papers published by Aeroballistics Division be tween 1959 and 1963. 

SECTION I: INTRODUCTION 

In support of the Army Ballistic Missile Agency's development work 
on the Jupiter missile and on satellites, the Fluid Dynamics Section of 
the Aeroballistics Laboratory was assigned in January 1958 to study the 
aerothermodynamics of ablation- type heat protection devices. With the 
exception of Reference 27, only experimental results and empirical re- 
lations for the effective heat of ablation were known at that time. For 
this reason, the purpose of the work to be done by the Fluid Dynamics 
Section was defined early in 1958 as follows: (1) to gain a better under- 

standing of the heat rejection mechanism in ablation- type shields, (2) to 
study the transient characteristics of this mechanism along flight tra- 
jectories, and (3) to determine the merits, disadvantages, and limitations 
of different types of ablative heat protection devices. Since initially 
only two and later four men were assigned to work on this problem, the 
ablation studies were restricted to the analysis of the transient heat 
and mass flows in homogeneous and noncharring materials which are mathe- 
matically tractable. Relations for heat and mass transfer at the ablating 
surface were supposed to be known from pertinent boundary layer solutions. 
The first step, therefore, consisted of the derivation of rigorous numeri- 
cal calculation methods which necessarily involve the solution of partial 
differential equations and thus demand computer evaluations. These calcu- 
lation methods were employed in 1959 and 1960 to study the performance of 
various glass and Teflon heat shields attached to several reentry bodies 
of revolution flying at zero angle of attack. The change of mission that 
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accompanied the transfer of this organization from ABMA to NASA on 
July 1, 1960, brought about a change of emphasis and left part of the 
initial program uncompleted. The continuation of the theoretical abla- 
tion studies by the Fluid Dynamics Section after July 1, 1960, was governed 
by the following purposes: (1) to augment' and finish work on some parti- 
cularly important phases of the initially outlined program of ablation 
studies, (2) to fill a request from the Jet Propulsion Laboratory for 
supporting work on the preliminary design of a Mars entry vehicle, (3) to 
fill a request from the Ames Research Center for rigorous numerical solu- 
tions pertaining to ablating glassy bodies, and (4) to fill a request 
from the Goddard Space Flight Center for supporting work on the aero- 
thermodynamics of tektite flight in order to contribute new aspects to 
this problem in space physics. 

SECTION II. CALCULATION METHODS 
A. Methods Derived and Employed 

The flows of heat and mass in objects ablating under conditions 
of hypersonic flight are governed by partial differential equations of 
the boundary layer type which represent the balances of mass, momentum, 
and heat. These equations are strongly nonlinear because of the tempera- 
ture dependency of the material properties. Table 1 lists the numerical 
solution methods which have been derived and applied for the case of 
homogeneous and noncharring materials.* The numerical methods (1) - (3) 
in the following are exact insofar as the numerical errors can be shown 
to tend to zero uniformly together with step sizes: 

(1) A transient difference method with two independent variables 
has been derived for the stagnation point of a melting and evaporating 
shield. This method is applicable at any cross section downstream of 

the stagnation point if flows of heat and mass parallel to the surface 
are negligible. 

(2) A transient difference method with three independent vari- 
ables has been derived for any cross section of a melting and evaporating 
shield of a body of revolution flying at zero angle of attack. 

These difference methods employ specially adapted procedures in the vicin- 
ity of the shield-air interface, where temperature gradients of the order 
of 1000 °K/mm are encountered during ablation of glass shields under condi- 
tions of hypersonic entry. The calculation of melt flow in method (2) is 
sensitive to discontinuities of the mathematical expressions for the gra- 
dients p'(x) or T^Cx) of pressure and interface shear stress, respectively, 
at the transition station xo of body segments possessing different geo- 
metric configurations, e.g., spherical nose and conical frustum. The 
following calculation method, with time as the only independent variable, 
is rigorously valid if the surface does not ablate and if heat conduction 
is restricted to the surface normal. 
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(3) The surface temperature is obtained from solving an integro' 
differential equation, whose solutions are good approximations to results 
of methods (1) and (2) in case of nonmelting but evaporating surfaces. 

Ablation type materials respond to high heating rates with high 
rates of surface radiation and/or mass transfer both in the molten and 
the gaseous states. For this reason, the net heat transfer rate across 
the surface is very much smaller than the gross amount of heat incident 
at the surface, provided this gross rate is sufficiently high. If this 
is true, the surface temperature may be approximated with a high degree 
of accuracy by ignoring both the preceding heating and ablation history 
and the net heat transfer rate as compared to the individual fluxes con- 
tributing to this net rate. Therefore, the following approximate calcu- 
lation method has been derived and applied: 

(4) The surface temperature is obtained from solving a time- 
independent ordinary equation which represents "radiation and mass trans- 
fer equilibrium." 

The initial program on ablation studies included the derivation of 
a Karman-Pohlhausen integral solution to serve as an engineering calcu- 
lation method for surfaces which melt and evaporate under transient condi 
tions.' This plan has not been carried out because of the transfer, from 
ABMA to NASA on July 1, 1960. 

B. Comparison of Transient and Quasi-Steady Calculation Methods 

Table 2 presents a list of typical objects whose descending 
flight across earth's atmosphere has been analyzed by use of the calcula- 
tion methods of Section IIA. Data for the stagnation points of the IRBM 
and the satellite appearing in Table 2 have been calculated by use of the 
transient method (1) in Section IIA and by use of the quasi-steady method 
presented in Ref. 1. Figure la shows the flight speeds, V, of these 
vehicles as functions of both altitude and entry time. Fig. lb presents 
the surface temperatures, Ti, Fig. lc the ablation speed v OT , and Fig. Id 
the distance a between the surface and the station where the temperature 
is 1000°F at the end of the flight. Since the quasi-steady analysis 
rests on a radiation and ablation balance of the incident heat, this 
method must fail before ablation begins. This is confirmed by the three 
performance graphs, Figures lb, lc, and Id. As a check on the accuracy 
of the calculation methods involved, one notices their close agreement 
during the ablation period. 

SECTION III. GENERAL CONCLUSIONS OF ABLATION STUDIES 

A. Review of the Problem 


The derived calculation methods have been employed to study the 
performance of shields made of glasses or of Teflon in case of selected 
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entry flights which cover the total range of missile, satellite, and space 
vehicle trajectories in the earth's atmosphere. In addition, the ablation 
process of stony meteorites has been studied numerically to supplement the 
field of small and moderate ablative mass losses of entry vehicles; this 
extension more clearly establishes trends in the variation of vehicle and 
trajectory parameters. 

The governing criterion for the heat protection of any vehicle tra- 
versing an atmosphere at hypersonic speed is the ratio of heat shield 
weight to the total weight of the vehicle. At minimum weight, a heat pro- 
tection device should keep both the temperature and the heat transfer at 
its inner surface below specified bounds. Therefore, as much of the 
incident aerodynamic and radiative heat transfer as possible should be 
returned into the air flow adjacent to the shield with the smallest 
possible ablative mass loss. 

If any object traverses an atmosphere in hypersonic, descending flight, 
the rates of both the aerodynamic and radiative heat fluxes incident at 
the surface are significant only when the flight speed is sufficiently 
high. Table 2 presents the altitude ranges in which aerodynamic heat 
transfer and ablation are encountered by typical objects traversing the 
earth's atmosphere. Because of the relatively small flight speed reached 
at the end of the heating period, heating and ablation occur only during 
a small fraction of the total flight time in the atmosphere. With the excep- . 
tion of cases like' short-range missiles, however, the magnitude of the peak 
heat transfer rate does not permit the utilization of heat shields only as 
heat sinks. 

Table 3 presents a list of available heat protection materials. 

The calculation methods listed in Section II have been applied to glasses 
as typical high-temperature ablators and to Teflon, a typical low-temperature 
ablator. In the following three sections, several particularly important 
aspects of ablation type heat protection shields are discussed. In section 
II IB , the merits and disadvantages of high- and low-temperature ablators 
are treated. Ablation by mass transfer and ablation by melt flow are com- 
pared in section IIIC. Section IIIB deals with transparent and nontrans- 
parent materials. 

B. High Versus Low Ablation Temperatures 

While available ablation type materials do not possess well 
defined temperatures at which the pertinent changes of state occur, tem- 
perature ranges always exist in which the ablation phenomena take place. 

If melt flow and evaporation of the heat shield's constituents occur only 
at high temperatures, the ablative mass loss is relatively small. For 
these so-called high-temperature ablators, radiation from the shield into 
the air may provide a significant decrease of the shield's heat absorption. 
This type of heat transfer back into the air is desirable, since it does 
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not involve any loss of shield material. High surface temperatures 
evidently require sufficiently small values of the thermal diffusivity 
and relatively short flight durations in order to avoid thick, and thus 
heavy, shields. If thermal stresses are disregarded, glass shields may 
survive heat pulses of extremely high magnitude and sufficiently brief 
duration. Glass-protected surfaces, therefore, can cope efficiently with 
the intense and short-lasting heat pulses characteristic of steep high- 
speed entry flights. Since, in general, the terminal flight period after 
the end of the ablation pulse is short on these trajectories, the conduc- 
tion of previously absorbed heat across the shield remains within toler- 
able bounds during this period. 

The performance of .a low-temperature ablator is characterized by 
approximately equal progress rates of the ablating surface and the thermal 
penetration of the shield. While the heated layer thus is thin, the abla- 
tive mass losses are relatively large and radiation emitted from the shield 
is insignificant in general. Since radiative heating incident at the sur- 
face in general is not absorbed by the gas released from an evaporating 
or sublimating surface, a low-temperature ablator is likely to fail 
entirely if exposed to radiative heating; i.e., catastrophic ablation 
rates may be expected then. (A high-temperature ablator, though, can 
cope efficiently with both aerodynamic and radiative heating modes provided 
the material is opaque to radiation.) These disadvantages are compensated 
for by the fact that some low-temperature ablators are considerably better 
thermal insulators than available high-temperature ablators; e. g. , Teflon, 
possesses a thermal diffusivity which is about one order of magnitude 
smaller than the one of glasses. Teflon shields, therefore, can cope with 
long-lasting aerodynamic heat pulses of low magnitude which, however, 
exceed the limit tolerable ofor nonablating radiation-cooled shields. 

This type of heating is experienced on nearly horizontal entry trajec- 
tories with small deceleration rates. 

C. Ablation by Mass Transfer Versus Ablation by Melt Flow 

The heat of evaporation or sublimation absorbed at the surface 
of a shield under conditions of aerodynamic, i.e., convective and/or 
radiative heating, constitutes a continuously distributed heat sink at 
tl^e shield-air interface. In addition, mass transfer reduces the con- 
vective heat transfer' in case of a laminar air boundary layer flow by 
influencing the velocity and temperature profiles. This beneficial 
effect is considerably smaller in case of turbulent boundary layer flow. 
Also, vapor ei, ^tted by an ablating surface in general is very nearly 
transparent to the radiation of the gas cap behind the shock wave. . 

The flow of melt in a layer adjacent to the surface of an 
ablating shield transports heat absorbed by the shield. The melt flow 
in general proceeds in the direction of the air flow adjacent to the 
shield and solidifies in regions with small heat transfer downstream 
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of the stagnation point. Instabilities of the melt flow, particularly 
where the deceleration of the entry vehicle causes a local reversal of 
the melt flow's direction, may cause a spraying of melt. While mass 
transfer decreases the aerodynamic heat transfer across the laminar 
boundary layer, melt flow in general removes hot melt from the vicinities 
of maximum heat input, and thus causes a more uniform distribution of the 
heat penetrating the shield. 

Some detrimental effects of mass transfer on calculated heat shield 
performance have been shown in References 12 and 13 for a quartz shield 
mounted on a space vehicle traversing. the earth's atmosphere under the 
conditions listed in Table 2. If it were possible to suppress evaporation 
at the stagnation point, the necessary shield thickness could be decreased 
from 76 to 64 mm at that point. This is explained by the calculated result 
that the nonevaporating opaque shield reaches a«maximum surface temperature 
of 4623 °K and removes 95 percent of the incident convective heat at the 
stagnation point by surface radiation, whereas the evaporating shield reaches 
a maximum temperature of 3101 °K and radiates only 28 percent back into the 
air . 


The evaporation of the quartz shield for this space vehicle reentry 
is very sensitive to increases in the aerodynamic heat transfer. When 
the heat transfer coefficient at the stagnation point of the evaporating 
shield is increased by a constant factor 2 during the entire flight, the 
total ablated thickness rises by a factor of 2.4, while the necessary 
thickness of the shield increases by a factor 1.65. When the heat trans- 
fer coefficient is increased by the same factor 2 at a distance of 4 
meters downstream from the stagnation point, where ablation is negligible, 
the necessary thickness of the shield rises by only 5 percent. If the 
heat transfer coefficient is increased by the factor 2 at the stagnation 
point of the nonevaporating quartz shield, the necessary thickness increases 
only by 8 percent. £ 

It may be concluded that the large mass transfer on the nose of this 
reentry vehicle keeps the surface temperature at such a comparatively low 
level that radiative cooling cannot govern the cooling balance. Decidedly 
contrary to mass transfer cooling, this radiative cooling possesses the 
desirable insensitivity to increases in the heat transfer incident at the 
surface. 


D. Transparent Versus Nontransparent Materials 

The effects of radiative heat exchange, in addition to the trans- 
portation of heat by convection and conduction in the shield, have been 
studied at the shield's stagnation point axis of symmetry by (1) consider- 
ing convective heating only, (2) ignoring reflection and scattering of 
radiation, and (3) neglecting any net radiative flux across a cylinder 
surrounding the stagnation point axis. Because of the last assumption, 
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calculation method 1 in Section II is applicable after adding pertinent 
terms to the energy balance, so that this differential equation then 
becomes an integro-dif ferential equation. The modified calculation method 
1 has been employed to study the influence of internal radiation in case 
of a space ship, whose flight data are listed in Table 2, and in case of 
cylindrical glass rods under axial heating and ablation conditions. The 
results of this study are presented in Table 4 and Figure 2. The experi- 
mental data for the Pyrex rod in Figure 2 have been obtained by courtesy 
of Dr. Dean R. Chapman of the Ames Research Center, NASA. Unfortunately, 
the absorption coefficient CL of the test specimen is not known. Calcula- 
tions were therefore carried out for the absorption coefficient of clear 
window glass, CC = 100, and the one of opaque material, a = co. The com- 
parison of experimental and calculated ablation data reveals satisfactory 
agreement, as is seen in Figure 2. 

The analysis of the calculated cases yields the following conclusions 
pertaining to changes of a, the only material parameter affecting radiant 
heat exchange: 

(1) In general, both the surface temperature and the surface 
temperature gradient decrease together with CL, causing lower temperatures 
in a layer adjacent to the surface and higher temperatures elsewhere. 

(2) Correspondingly, when Ot is changed, the rate of melt flow 
may increase or decrease, whereas the evaporation rate changes together 
with CL. 

(3) The time-dependent length c(t) is defined as the distance 
between the surface and that station in a semi-infinite shield where any 
prescribed temperature is reached. This parameter c(t) can be used to 
determine the necessary thickness of the shield. Item (1) of this list 
shows that c(t) may increase or decrease upon any change of a. If the 
heating pulse stops a sufficient time before the end of the flight 

(t = t £ ) , the parameter b = c(t e ) increases as CL is reduced. 

This list is confirmed by the results presented in Table 4. It is seen 
that both the residual thickness b of the shield and its total necessary 
thickness (s + b) at the stagnation point of the space ship increase as 
CL decreases. However, (s 4- b) decreases together with CL in case of the 
rods made of Pyrex glass and quartz, where b and s pertain to the time 
instant the heating pulse terminates. 

The so-called effective, apparent, or fictitious thermal conductivity 
has been employed frequently in the past to calculate superimposed radia- 
tive and conductive heat fluxes by use of a substituting heat conduction 
mechanism. Mathematical solutions and experiments express the effective 
thermal conductivity under stationary conditions as dependent on the 
geometry of the problem and as proportional to the third power of the 
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absolute temperature. Because of the high temperatures encountered and 
the transient character of heat shield performance under entry conditions, 
use of the effective thermal conductivity gives grossly deviating results 
as compared to studies properly accounting for conductive and radiative 
transport of heat. 

E. The Influence of Individual Glass Properties on Shield Weight 

The necessary thickness of the glass needed on the shield to 
protect the vehicle is the sum of the total ablation s and the residual 
thickness b needed as a heat sink between the final surface and the 
supposedly thermally insulated inner edge. This thickness b is determined 
by the restriction that the thermally insulated inner edge never reaches a 
temperature larger than 460 °K at any time during the reentry. The neces- 
sary weight y(s + b) is the weight per unit area of the required glass 
shield, where 7 is the specific weight of the glass. 

The effect of the material properties of glass on the necessary 
weight was studied for many glass shields, each with a different set of 
assumed properties, in case of the IRBM, ICBM, and the satellite listed 
in Table 2. Figures 3 and 4 present some results of this study. 

The thermal properties of the glass shields were varied as follows: 

-4 -4 

(1) thermal conductivity k, 0.61 x 10 < k < 9.713 x 10 

(kcal/m°K sec), 

(2) specific heat Cp, 0.10 < Cp < 0.58 (kcal/kg°K), 

(3) emissivity constant e, 0.10 < e <0.8 (-), 

(4) viscosity level of Pyrex glass ji , 0.10 < \i < 100 (-)> and 

(5) vapor pressure level of Pyrex glass A, 0 < A < 1000 (-). 

A and p. are constants that the vapor pressure and viscosity functions 
p v (T) and p(T) , respectively, of Pyrex glass were multiplied by in order 
to study the effect of these material parameters. 

Experimental results for opaque Pyrex glass show thatk./7 c p the 
relevant combination of k, 7 , and Cp, slightly decreases with increasing 
temperatures; e.g. , T = 473 °K yields k/yc = 7.9 x 10" 7 m s /sec and T = 1073 °K 
yields k/yCp = 7.3 x JL0 " 7 m 2 /sec. The reiults presented in Figures 3 and 4 
were calcuated by assuming k, 7 , and Cp as constants. In order to check 
the effect of this assumption on k, 7 , and Cp, the necessary weight of a 
quartz shield for the space ship reentry listed in Table 2 was also calcu- 
lated by using curve fits for the material properties k = k(T) and c p = c p (T) 
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The results have been compared with a case where mean constant values for 
the two properties k and c p were employed. Both cases yielded approximate- 
ly the same value of the necessary weight. 

Figure 3 shows the total ablation s, the heat sink thickness b, and 
the necessary thickness (s + b) for the IRBM reentry as functions of k, 

Cp, and the viscosity level factor p* for Pyrex glass. It is evident 
from this graph that small values of k, large values of c p , and large 
values of p* are desirable properties of glass shields in order to mini- 
mize the necessary weight. Similar graphs for the ICBM and the satellite 
show the same trend. 

Figure 4 shows the effect of the material properties on the necessary 
weight at the stagnation point of the glass shields for three reentries 
listed in Table 2 as a function of the thermal diffusivity k/yc p . The 
shaded areas show for each reentry the range of necessary weights that 
result from varying the thermal properties between the previously mentioned 
limits. It is evident from this graph that k /yc' is £ he most influential 
material parameter on the necessary weight for all of the investigated 
reentries if k/yc has the desirable small values and/or ablation is insign 
4ficant. This is P verified by the calculated results which show for given 
material properties that the IRBM experiences the smallest and the ICBM 
the largest ablative mass loss. As k/yc p increases, the temperature pro- 
file normal to the surface spreads, and thus the molten layer becomes 
thicker. Correspohdingly , the influence of each individual material pro- 
perty k, c p , e, p, and p v on y(s + b) increases together with k/yCp. 

Numerical solutions for heat shield performance on ballistic entry 
trajectories show that large percentages of both the time-integrated sur- 
face radiation and the time -integrated mass transfer rates are concentrated 
in the short-lasting, high-temperature ranges. According to these calcu- 
lations, changes of' all the relevant material properties have only very 
little effect on the surface temperature in this range because of the 
strongly nonlinear temperature dependencies of the material properties and 
of the Stefan -Boltzmann radiation law. For a given thermal diffusivity, 
then, the temperature profile normal to the surface of an opaque shield 
depends only weakly on changes in e, p v , and p. 

The numerical solutions also indicate that any increase of ablative 
mass losses by melt flow or evaporation due to property changes are com- 
pensated for to a large extent by a smaller heat absorption in heat sink 
fashion and thus also a smaller thermal penetration of the shield. This 
explains, e.g. , why changes of the vapor pressure level A in the bounds 
0 < A < 10 3 influence the necessary shield thickness (s + b) only within 
the limits ± 7 percent, as solutions for the four reentry vehicles listed 
in Table 2 show. 

SECTION IV. CALCULATED RESULTS FOR INDIVIDUAL RE-ENTRY FLIGHTS 

A. Space Ship Traversing the Atmosphere of the Earth 

The performance of quartz and Teflon shields was studied for a 
space ship traversing the atmosphere of the earth (see Ref. 12). The 
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capsule is presented in Figure 5a, the entry parameters appear in Table 2. 
The capsule possesses a spherical nose of one meter diameter followed by 
a conical frustum of 3.69 meters in length with a largest diameter of 3.77 
meters. This capsule has a volume of 18 m 3 and weighs 8640 kg. The 
ballistic factor is 500 lb/ft 2 under the assumption of a mean drag coeffi- 
cient cp of 0.32. The reentry point was placed at the altitude of 120 km 
above sea level where both aerodynamic and radiative heating are still 
negligible. 

Calculations have been made at the stagnation point for both a compact 
quartz shield and one with 50 percent porosity, hereafter called foamed 
quartz shield. While the specific weights differ by a factor 2, the ther- 
mal dif fusivities k//Cp were assumed equal. Also, transportation of heat 
by radiation across the inner cavities was excluded. The results show 
that both shields experience approximately the same weight loss per unit 
■area at the stagnation point. In order not to exceed the prescribed tem- 
perature of 460°K alt the shield's insulated inner edge, the compact shield 
needs to have a thickness of 44 mm' while the foamed shield needs only 39 mm 
because the surface temperature distributions are approximately equal, and 
the smaller thermal conductivity of the foamed material causes larger sur- 
face temperature gradients. The necessary weight per unit area of the 
compact shield at the stagnation point therefore exceeds the one of the 
foamed shield by a factor 1-66, in order to minimize the change in the 
vehicle's shape, the ablating outer layer of the quartz shield should be 
compact, while the inner layer must be porous for small shield weight. 

The performance of a foamed quartz shield and a Teflon shield have 
been analyzed on the spherical and the conical parts of the capsule sur- 
face. The calculations indicate that 99 percent of the quartz's ablative 
mass loss is due to evaporation so that the calculation method No. 1 in 
Section IIA is applicable for the entire body. Figures 5b and 5o, 
present for the quartz shield the surface temperature and the vapor flow, 
respectively, as functions of time and distance x from the stagnation 
point. The subsidiary graph in Figure 5b shows the time when the assumed 
transition Reynolds number of 10 s is reached at station x tr . The second 
peak of heat transfer and surface temperature occurs immediately after 
transition, and results in negligible ablative mass of the quartz shield. 

On the conical frustum of the Teflon shield, however, 40 percent of the 
total mass loss takes place after transition has been reached. 

Figures 6a and 6b present the surface temperatures at the stagnation 
point of the foamed quartz shield and the Teflon shield, respectively, as 
functions of time t and distance y* normal to the surface. It is seen that 
Teflon ablates considerably faster than quartz at the stagnation point. 
Because of both a smaller thermal diffusivity and a more rapid rate of 
surface recession, the temperature profile in the Teflon shield is consider 
ably steeper than the one in the quartz shield. Accordingly, the Teflon 
shield absorbs significantly less heat than the quartz shield. 
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Figure 7 shows the necessary weight of the quartz shield and the 
Teflon shield as functions of x. Since the heat shield performance on 
the nose is governed by the ablative mass loss, the quartz shield weighs 
less than the Teflon shield on the nose. Because the heat absorption and 
thermal penetration are dominant on the frustum, the Teflon shield is 
lighter than the quartz shield there. 

The effects of the thermal radiation incident at the surface of the 
space ship's heat shield have been ignored up to. this point. Available 
pertinent data indicate that the gas radiation incident at the stagnation 
point amounts to as much as 65 percent of the aerodynamic heat transfer 
to a cold wall for the reentry flight under discussion. If the quartz 
shield is nontransmittent to thermal radiation, its necessary thickness 
at the stagnation point increases by only 16 percent; this is explained 
by the fact that the maximum surface temperature at the stagnation point 
rises only from 3101 °K to 3161 °K when the gas radiation is added to the 
aerodynamic heat transfer at the stagnation point of the quartz shield. 

B. Capsule Traversing the Atmosphere of Mars 

The unmanned Mars probe consists of a spherical nose and a 
conical afterbody. The maximum diameter is 0.814 m and the overall 
length 1.28 m. The entry parameters at the altitude of 300 km above 
Mars' surface (entry time zero) are listed in Table 2, The flight velo- 
city decreases from 6.5 to 1 km/sec in the period 40 to 80 sec, followed 
by a long period of low-speed flight to impact at 210 sec. This gives a 
short pulse of intense aerodynamic heating with a peak of 500 kcal/m 2 sec 
at the stagnation point during the period 30 to 80 sec, followed by a 
long period of negligible heating. Thus, a material for the heat shield 
has to be found which minimizes the sum of ablated thickness s and residual 
shield thickness b at impact time. 

Four materials were considered for the ablation type heat shield: 
quartz, fibergla.s, Teflon, and phenolic nylon. For quartz and fiberglas 
the calculation method no. 2 in section II A has been used, where fiber- 
glas shield was approximated by a homogeneous model. Teflon and phenolic 
nylon decompose without forming a molten layer. A modification of calcu- 
lation method no. 2 was needed to account in an approximate way for the 
formation of a char layer on an ablating shield made of phenolic nylon. 
Experimental results for the effective heat of ablation and the ablation 
temperature were used for the phenolic nylon calculations. 

Even considering the most conservative of the available empirical 
expressions, phenolic nylon was found to make the lightest shield. The 
main advantages of phenolic nylon are its low thermal diffusivity, the 
high emissivity of the char formed by decomposition, and the fact that 
the decomposition forms gases of low molecular weight, which are: more 
effective in heat blocking than gases of high molecular weight. Quartz 



12 


was the least effective shield material because it did not have enough 
vaporization for heat blocking at the temperatures reached, and its thermal 
diffusivity was so high it allowed a greater thermal penetration than the 
other materials. 

The structure under the shield was supposed to act as a heat sink, 
thus helping the shield dissipate the absorbed heat. Two structure materials 
were considered, aluminum and fiberglas. When the shield is thick, the heat 
flux penetrating to the structure is small and distributed over a long time 
interval. This makes the structure's heat capacity more important than its 
thermal conductivity. On the conical section of the capsule, where only a 
thin shield is needed, the heat flux to the structure is in a short pulse. 

A high thermal conductivity of the structure is mandatory on the cone to 
dissipate the heat quickly. From these thermal considerations aluminum 
was found to make the lighter structure. 

SECTION V. FLIGHT OF NATURAL QUARTZ OBJECTS IN THE EARTH'S ATMOSPHERE 

A. Tektite Flight 

Tektites are small glassy bodies which are found in Indochina, 
the Philippines, Indonesia, Australia, Texas, Georgia, Ghana, Bohemia, 
and Moravia. All tektites have a family-like chemical composition and 
are unrelated to the local geological formations; therefore, they must 
have been hurled up somewhere by violent natural events and carried in 
flight into the strewn fields. All the australites and some javaites 
belong to one of a few well-defined geometric shapes, which can be 
explained by ablation in hypersonic flight in the earth's atmosphere. 

It has not been possible yet to prove whether tektite flight started at 
the eartHs surface or at some extraterrestrial point of departure. 

Chapman has presented conclusive experimental proof in Ref. 25 that glass 
spheres, when placed in an electric-arc jet tunnel and exposed to heating 
rates of the order of those experienced in hypersonic flight, become 
strikingly similar to the button-type australites (see Fig. 8). Chapman 
also has shown in Ref. 25 that an oscillation about any axis other than 
the flight axis is rapidly damped by the continuously increasing magni- 
tude of the pressure distribution in descending flight, whereas this 
damping effect is absent in ascending flight. 

Since the button-type australites then must have entered the earth's 
atmosphere as cold glassy spheres, ablation an<^ trajectory analysis are 
applicable. Because the temperature level in tektites may take any values — 
whereas the temperature and the temperature gradient at the inner edge of 
the heat shield of a vehicle are subjected to upper bounds — the conclusions 
from the tektite study differ somewhat from the ones of vehicle studies. 

The results for the flight of the button-type australites to be presented 
in the following have been obtained under the assumption of a simplified 
geometrical model whose shape is determined by the stagnation point ablation 
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only. Experimental results show that radiative heating of objects as 
small as tektites is negligible in hypersonic flight. Fig. 9 presents, 
as functions of entry speed and entry angle 0^, curves of constant 
ratio of total heat blocked by mass transfer over total aerodynamic heat 
flux, both at the stagnation point. It is seen' that more than 90 percent 
of the aerodynamic heating rate is blocked at high entry speeds. Calcu- 
lations show that less than 20 percent of the total kinetic energy converted 
into heat enters the boundary layer surrounding the object. These consider- 
ations indicate why in Figure 10, which presents curves of constant relative 
mass loss as functions of entry speed and entry angle, spheres of only 11 
grams initial weight may enter with as much as 26 km/sec in the vertical 
direction and still lose only 50 percent of their initial mass. (Both 
convective and radiative heating rates at entry speeds i > 15 km/sec have been 
obtained by extrapolating correlation formulas valid for lower speeds.) 

Table 5, which has been compiled by analyzing all the calculated 
results for tektite flight, explains the dependence of relevant parameters 
on entry speed Vj_, entry angle 0 i , initial mass m(0) , and the level of the 
vapor pressure of the material. This dependence is predominantly deter- 
mined by the aerodynamic heating pulse and by the shield's mass transfer 
effect* A smaller percentage of the kinetic energy converted into heat 
reaches the surface as the heat pulse is shifted to lower altitudes, where 
a higher portion of the generated heat remains in the denser air. The 
case is clearest for the change of the entry angle 0^, since this does 
not affect the total kinetic energy converted into heat. The following 
conclusions follow from Table 5: 

(1) The relative mass loss increases as the initial mass m(0) 
decreases; i.e. , the chance of survival drops together with the initial 
size of small objects. 

(2) The relative mass loss increases only moderately together 
with the entry speed because the heat pulse is shifted to a higher alti- 
tude so that the mass transfer effect increases more rapidly than the 
surface temperature does. This explains why small objects can traverse 
the atmosphere at high entry speeds. 

(3) As the trajectory becomes steeper (0^ increases), the 
relative mass loss drops significantly since the time -integrated aero- 
dynamic heating pulse, calculated as if the surface were cold, decreases 
strongly; from this it is seen that the chance of survival of small bodies 
is better on a steep than on a shallow--i.e. , long-lasting--trajectory . 

(4) As the level of the vapor pressure increases, a strong decrease 
of both the surface temperature and the relative mass loss takes place, 
while the mass transfer effect rises significantly. However, this decrease 
in relative mass loss is not true, in case of a reentering nose cone where 
the heated quartz layer is relatively small. 
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B. Meteor Flight 

Since the thickness of the hot gas cap between the shock front and 
the surface of an object traversing the atmosphere in hypersonic flight 
increases together with the body's cross section, meteors are predominantly 
heated by equilibrium or nonequilibrium radiation from the hot gas cap, 
whereas bodies as small as tektites essentially experience convective 
heating only. Detailed calculations show that the extremely high radia- 
tion rates of as much as 10 4 - 10 s kcal/m 2 sec, which are incident at 
the surface of a spherical meteor- of 10 s tons, cause the surface to boil. 

If the material is opaque to thermal radiation, the heated layer beneath 
the rapidly receding surface has a thickness of considerably less than 
1 mm. The calculations show, however, that a layer of several centimeter 
thickness reaches temperatures of 3000 to 5000°K if the material is suffi- 
ciently transparent to thermal radiation. Vapor nuclei are then formed 
at nonhomogeneous spots in the glassy melt. Unless strong temperature 
gradients normal to the surface exist in the melt, there is no force present 
which drives the bubbles towards the surface. The outer portions of the 
radiatively heated layer therefore consist of a spongy glass melt, which 
offers little resistance to being removed from the meteor by the air 
forces acting in hypersonic flight. It may be concluded that the intensity 
of the gas radiation rates, which is comparable only to the ones liberated 
by nuclear explosions, causes any heavy object traversing the atmosphere 
in hypersonic flight to experience large mass loss rates regardless of 
the material's transparency to thermal radiation. 



TABLE I 


LIST OF CALCULATION METHODS 


Number 

of 

\ Method 

Number of 
Independent 
Variables 

Number of 
Space 
Variables 

Transient 

or 

Stationary 

Method Furnishes 
Numerically 
Exact Solution Of 

Characteristics of 
Numerical Solution 

1 

' 

2 

■ 

Transient 

3 Partial Differen- 
tial Equations 

Modified Difference 
Method 

2 

3 

2 

Transient 

3 Partial Differen- 
tial Equations 

Modified Difference 
Method 

3 

1 

0 

Transient 

Integro -Differential 
Equation for Surface 
Temperature 

Step Method Pro- 
ceeding in Time 
Direction 

4 

0 

0 

Stationary 

Ordinary Equation 
For Surface 
Temperature 

. 

Solution is 
Determined by 
Graphical or 
numerical 
methods 


Solution Yields 

General Remarks 

Temperature and Melt 
Flow Profiles, Sur- 
face Evaporation or 
Sublimation Rates 

Applicable 

(a) in the vicinity of a 
stagnation point with or 
without radiant heat 
exchange 

(b) everywhere if the 
flows of heat and mass 
parallel to the surface 
are negligible and if 
radiant heating may be 

. ignored. 

Temperature and Melt 
Flow Profiles, Sur- 
face Evaporation or 
Sublimation Rates 

Applicable to any melting 
and evaporating shield of 
a body of revolution fly- 
ing with zero angle of 
attack if the shield is non- 
transmittent to radiation. 

Surface Temperature 
and Evaporation or 
Sublimation Rates 
at Surface 

(a) Method furnishes an 
exact solution for nonablat; 
ing surfaces and one- 
dimensional heat conduction 
without radiant heat 
exchange in the shield. 

(b) Method yields good 
approximation up to the end 
of the ablation period if 
melt flow and radiant heat 
exchange in the shield are 
absent 

Surface Temperature 
and Evaporation or 
Sublimation Rates 
at Surface 

Method is valid if the net 
heat flux across the sur-^ 
face is very much smaller 
than the individual contri- 
butions to the net heat 
flux and in the absence of 
melt flow and radiant heat 
exchange . 




























TABLE II 


ALTITUDE RANGES OF AERODYNAMIC HEATING AND ABLATION PULSES 


Entering 
Object as 
a Typical 

Entry Angle 
Relative to 
Earth's 

Entry 

Velocity 

(km/sec) 

1 

'At the Stagnation Point 




Aerodynamic Heating Values Pertain to Wall at 

300 °K 

If Radiation of 

Shock Layer is 
-Ignored 

Represen- 
tative of 

(Mar 1 s ) 
Horizon 


Aerodynamic 
Heating Begins 
to exceed 100 
Kcal/m 2 sec at 
Altitude (Km) 

Aerodynamic 
Heating 
Reaches its 
Maximum at 
the Altitude 
(Km) 

Aerodynamic 

Heating 

Becomes 

Smaller 

Than 100 Real/ 
m 2 sec at 
Altitude (Km) 

Maximum 
Value of 
Aerodyna- 
mic Heat- 
ing (Kcal/ 
m 2 sec) 

Ablation 
of Quartz 
begins at 
Altitude 
(Km) 

Ablation 
of Quartz 
Shield 
Ends at 
Altitude 
(Km) 

IRBM 

35° 

4.5 

67.5 

25 

12 

840 

36 

10 

ICBM 

24° 

7.226 

86.5 

20 

5.9 

6500 

71 

6 

Satellite 

2° 

7.9 

90 

55 

37 

625 

76 

36 

Space Ship 

5° 

11 

96.6 

47 

23 

2800 

80 

22 

Tektite 

(Meteorite 

. 

6° 

7 

90.0 

68 

, 56 

1000 

75 

51 

Tektite 

(Meteorite 

90° 

30 

106 

61 

37 

248,895 

92 

37 

Mars 

Vehicle . 

40 

6.4 

145 

62 

40 

500 

110 

40 




TABLE III 


CLASSES OF ABLATORS 



HT = High temperature ablator 
LT = Low temperature ablator 








TABLE IV 


EFFECTS OF INTERNAL RADIATION AT STAGNATION POINT AXIS 


Space Capsule with Quartz Shield Treated in Section IIIA 



Absorption Coefficient 

1/m 

100 

200 

00 

Maximum Surface Temperature at Stagnation Point 

°K 

2974.5 

3023.0 

3100.0 

Evaporation Rate at Time of Maximum Surface Temperature 

mm/sec 

.0518 

.0684 

.1135 

Distance c Between Surface and Station 460°K at Time of 
Maximum Surface Temperature 

mm 

39.85 

45.0 

2.013 

Time-Integrated Ablation Rate 

mm 

5.46 

9.48 

19.117 

Final Thickness b of shield 

mm 

81.01 

63.0 

43.86 

Necessary Thickness (s + b) of Shield 

mm 

86.47 

72.48 



62.98 



Results for Axial 
Heating of Pyrex 
Cylinder 

Results for Axial 
Heating of Quartz 
Cylinder 

Absorption Coefficient 

1/m 

100 

00 

100 

00 

Time Instant = Total Test Time 

sec 

10 

— 

10 

10 

10 

Surface Temperature at t = 10 sec 

°K 

2068 

2050 

2589 

'2622 

Distance b Between Surface and Station 460 °K at t = 10 sec 

mm 

3.25 

3.33 

5.5 

6/5 

Time Integrated Calculated Ablation Rate s 

mm 

3.92 

4.55 

00 

.55 

•k 

Time Integrated Experimental Ablation Rate s x 

mm 

3.98 


.39 


Distance (s + b) of Station 460°K from Original Surface 

mm 

7.17 

7.88 

5.98 

7.05 





TABLE V 

GENERAL CONCLUSIONS DRAWN FROM THE ANALYSIS OF TEKTITE FLIGHT 
The heat transfer parameter listed and the temperature pertain to the stagnation point. 


Parameter 

Change due to Increase 
of Entry Speed 

Change due to Increase 
of Entry Angle 

Change due to Increase 
of Mass at Entry Time 

Change due to Increase; 
of Vapor Pressure 
Level 

Relative Mass Loss 

moderate increase 

strong decrease 

moderate decrease 

strong decrease 

Aerodynamic Heat 
Transfer times 
Crossectional Area 
over kinetic energy 
converted into heat 

small increase 

small decrease 


very small increase 

Altitude where 
Maximum of Aero- 
dynamic Heat Pulse 
occurs 

small increase 

moderate decrease 

small decrease 

no change 

Duration of Aero- 
dynamic heat pulse 


strong decrease 

moderate increase 

no change 

Maximum surface 
temDerature 

small increase 

moderate increase 

moderate increase 

strong decrease 

Maximum Aerodynamic 
heating rate 

strong increase 

strong increase 

moderate increase 

very small increase 

Ratio of heat radia- 
ted from surface to 
Aerodvnamic neatine 

moderate decrease 

moderate decrease 

nearly constant 

strong decrease 

Ratio of heat blocked 
by mass transfer to 
aerodynamic heating 

moderate increase 

small increase 

very small increase 

strong increase 

Ratio of Evaporation 
to total ablation 

moderate increase 

small increase 


strong increase 

















































Distance Ablated (mm) 


0 Experimental Data from Ames 
O Calculated with Internal Radiation («= 100) 
^Calculated without Internal Radiation 




Time (sec) 


FIGURE 2: COMPARISON OF EXPERIMENTAL AND CALCULATED RECESSION OF SURFACE AT STAGNATION POINT 

OF PYREX CYLINDER 



PLANES REPRESENT 


to 

to 


/A* = 1.0 


A a* = 0.1 
O LL* =1.0 



FIGURE 3 MAXIMUM THERMAL PENETRATION b AND TOTAL ABLATION s AT 
THE STAGNATION POINT OF AN IRBM 

















FIGURE 5 A. SURFACE TEMPERATURE T[x,o,3 FOR A QUARTZ SHIELD AND 

TRANSITION POINT x tr [t]OF THE AIR BOUNDARY LAYER (SPACE SHIP) 
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X (m) 


FIGURE 7: NECESSARY WEIGHT OF GLASS AND TEFLON SHIELDS AS FUNCTIONS OF 

DISTANCE X FROM STAGNATION POINT FOR RE-ENTERING SPACE VEHICLE 
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FIGURE 8: COMPARISON OF THREE AUSTRALIAN BUTTON TEKTITES (AT RIGHT, NOW 

IN BRITISH MUSEUM) WITH THREE TEKTITE GLASS MODELS (AT LEFT) 
ABLATED BY AERODYNAMIC HEATING 
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Data pertains to initially hemispheric model with: 

Radius R(0) = 1.30 cm, i.e. W(0) = 11.04 grams, and vapor pressure p vl (T) ( ) 

Radius R(0) = 0. 65 cm, i.e. W(0) = 1.38 grams, and vapor pressure p vl (T) ( ) 

Radius R(0) = 1.30 cm, i.e. W(0) = 11.04 grams , and vapor pressure p V2 (T) ( ) 

q = (q -q )+vhV, where q - q is due to the diffusion of vapor 

H bi \ero aero v s aero \ero v 

across the boundary layer and yh v V s is the heat absorbed by the evaporation process. 




HEATING,/ q w d1 /f q aero dt, AT THE STAGNATION POINT 
oo 

(AUSTRALIAN BUTTON TEKTITE) 





Entry Velocity Vi(km/sec) 


50 


Data pertains to vapor pressure function p (T) and initially hemispheric 
model with mass m(0): 

Radius R(0) = 1.30 cm, i.e. W(0) = 11.04 grams ( ) 

Radius R(0) = 0 . 65 cm, i.e. W(0) = 1 . 38 grams (— ) 

m{tf) = final mass at impact time 



FIGURE 10: RELATIVE MASS LOSS, l-m(t f )/m(0) OF BUTTON TYPE 

AUSTRALITES 
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